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MAPPING THE STABILITY REGION OF THE 3:2 NEPTUNE-PLUTO RESONANCE 

II. F. Levison k S.A. Stern (Southwest Research Institute. San Antonio. TX) 

Pluto and Charon are most likely the remnants of a large number of objects that existed in the 
Oamis-Neptune region at early epochs of the solar system [1]. Numerical integrations have shown that, 
in general, such objects were ejected from the planetary region on timescales of % 10' years (cf. [2]) after 
Neptune and Uranus reached their current masses. It is thought that the Pluto-Charon system has survived 
to current times without being dynamically removed in this way because it is trapped in a set of secular 
and mean motion resonances with Neptune (cf. [3]). The best-known Pluto-Neptune orbit coupling is the 
3:2 mean motion resonance discovered almost 30 years ago by C. Cohen and E. Hubbard [4]. These workers 
showed that the resonance angle, 8 = 3A p — 2X>\ — Sjp. where 2'p is the longitude of perihelion of the 
Pluto-Charon system, and A,v and A p are the mean longitude of Neptune and Pluto-Charon respectively, 
librates about 180° with an amplitude, As>, of 76°. 

We report here on a numerical simulation project to map out the stability region of the 3:2 resonance. 

The results of these simulations are important to understanding whether Pluto’s long-term heliocentric 
stability requires only the 3:2 resonance, or whether it instead requires one or more of the other Pluto-Neptune 
resonances. 

Our study also has another important application. By investigating stability timescales as a function 
of orbital elements, we gam insight into the fraction of orbital phase space which the stable 3:2 resonance 
occupies. This fraction is directly related the probability that the Pluto-Charon system (and possibly other 
small bodies) could have been captured into this resonance. 

In the simulations we report here, we employed a numerical integrator to evolve t he orbits of massless 
test particles in and around the 3:2 resonance under the direct gravitational influence of Uranus and Neptune. 

The test particles were followed using the simplectic integration scheme described in [5]. The motions of 
Uranus and Neptune were calculated analytically from a synthetic secular perturbation theory using the 
four dominant perihelion precession frequencies discussed in [0]. Thus, although Jupiter and Saturn were 
not explicitly included in the simulation, their effect on the orbits of Uranus and Neptune are included. This 
approach has been shown to give good statistical agreement with techniques that include direct integration 
of the planets for particles in this region [7]. The integrations lasted 2 billion years. The particles were 
followed until they suffered a close approach to either Neptune or Uranus, at which time they were assumed 
to be removed from the resonance. 

Our goal is to determine whether the orbit of a test particle is stable over a significant fraction of the 
age of the solar system as a function of As, and its initial inclination and argument of perihelion. /* 0 and u>n 
respectively. VVe plan to study orbits with As < 180°, / 0 < 180°. 0 < u; 0 < 90°, and an initial eccentricity 
equal to the current eccentricity of Pluto, 0.248. 

We present here a preliminary report on this study. Over 160 test particles have been integrated in 
order to sample important aspects of the orbital element space. As of this time we have studied the regions 
of As < 180°, io < 35°, with u either 0 or 90°. 

First we present the results for our most well studied case, particles with iq = 17.2°, which is the 
inclination of Pluto. Figure 1 shows the removal times for these particles as a function of their libration 
amplitude. As - The open and filled circles represent objects with cjo = 0 and 90° respectively. If a symbol is 
on the dotted line then the particle survived for the length of the integration (2 billion years). The largest 
. 4,5 that survived is approximately 80°. Interestingly, the Pluto-Charon system has 4 , 5 = 76°. 

An examination of the individual orbits shows that the argument of perihelion circulates rather than 
librates. Since uj of our particles circulate, there should be no fundamental difference between the uj = 0 and 
uj — 90° runs, and this is what is observed. The argument of perihelion of the Pluto-Charon system librates 
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around 00° so that it i.s above the orbital plant' of Neptune when it is al perihelion ( [$] ) . It has been argued 
that this aids m the stability of its orbit. Our results show that this may not ho necessary. Indeed, we lind 
that objects in the 3:2 resonance with Neptune that have the same A.\ and inclination as the Pluto-Charon 
system seem to be stable for at least 2 billion years, even if they are not protected by the libration in w. 

Preliminary results for objects with different inclinations are as follows. We find that the maximum 
libration amplitude that, is stable. A^ Miax is not a strong function of either i {) or wq. For most values of / (J 
and wo, . U max ~ 90°. The only exception appears to be for a small region near u; t) = 90° for /'n ^ 20°. For 
this region all orbits appear to be short lived, having lifetimes less than 10 8 years. The explanation for this 
is being investigated. 


00 

o 

hJ 





Figure 1 - The removal times of massless test particles in the 3:2 resonance with Neptune. There initial 
inclination is 17.2°, the initial argument of perihelion is either 0 (open circles) or 90° (filled circles). 
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